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We have searched for evidence of physics beyond the standard model in events that include an 
energetic photon and an energetic 6-quark jet, produced in 85 pb~^ of pp collisions at 1.8 TeV at 
the Tevatron Collider at Fermilab. This signature, containing at least one gauge boson and a third- 
generation quark, could arise in the production and decay of a pair of new particles, such as those 
predicted by Supersymmetry, leading to a production rate exceeding standard model predictions. 
We also search these events for anomalous production of missing transverse energy, additional jets 
and leptons (e, fi and r), and additional 6-quarks. We find no evidence for any anomalous production 
of 76 or 76 + X events. We present limits on two supersymmetric models: a model where the photon 
is produced in the decay X2 ~^ 7Xij ^-nd a model where the photon is produced in the neutralino 
decay into the Gravitino LSP, Xi ~^ jG. We also present our limits in a model-independent form 
and test methods of applying model-independent limits. 

PACS number(s): 13.85.Rm, 13.85.Qk, 14.80. Ly 



I. INTRODUCTION 



As the world's highest-energy accelerator, the Tevatron Collider provides a unique opportunity to search for evidence 
of physics beyond the standard model. There arc many possible additions to the standard model, such as extra 
spatial dimensions, additional quark generations, additional gauge bosons, quark and lepton substructure, weak-scale 
gravitational effects, new strong forces, and/or supersymmetry, which may be accessible at the TeV mass scale. In 
addition, the source of electro-weak symmetry breaking, also below this mass scale, could well be more complicated 
than the standard model Higgs mechanism. 

New physics processes are expected to involve the production of heavy particles, which can decay into standard 
model constituents (quarks, gluons, and electroweak bosons) which in turn decay to hadrons and leptons. Due to 
the large mass of the new parent particles, the decay products will be observed with large momentum transverse to 
the beam (pt), where the rate for standard model particle production is suppressed. In addition, in many models 
these hypothetical particles have large branching ratios into photons, leptons, heavy quarks or neutral non- interacting 
particles, which are relatively rare at large values of pt in ordinary proton-antiproton collisions. 

In this paper we present a broad search for phenomena beyond those expected in the standard model by measuring 
the production rate of events containing at least one gauge boson, in this case the photon, and a third-generation quark, 
the b-quark, both with and without additional characteristics such as missing transverse energy (^J. Accompanying 
searches are made within these samples for anomalous production of jets, leptons, and additional b-quarks, which are 
predicted in models of new physics. In addition, the signature of one gauge-boson plus a third-generation quark is 
rare in the standard model, and thus provides an excellent channel in which to search for new production mechanisms. 

The initial motivation of this analysis was a search for the stop squark (t) stemming from the unusual ee"f"f event 
observed at the Collider Detector at Fermilab (CDF) A model was proposed |^ that produces the photon from 
the radiative decay of the X2 neutraHno, selected to be the photino, into the x?, selected to be the orthogonal state 
of purely higgsino, and a photon. The production of a chargino-neutralino pair, xt^2^ could produce the jbfJ^ final 
state via the decay chain 

X+X§-(?fo)(7X?)-(6cx?)(7X?). (1) 

This model, however, represents only a small part of the available parameter space for models of new physics. 
Technicolor models, supersymmetric models in which supersymmetry is broken by gauge interactions, models of new 
heavy quarks, and models of compositness predicting an excited b quark which decays to 76, for example, would also 
create this signature. We have consequently generalized the search, emphasizing the signature (76 or jb^^) rather 
than this specific model. We present generalized, model-independent limits. Ideally, these generic limits could be 
applied to actual models of new physics to provide the information on whether models are excluded or allowed by the 
data. Other procedures for signature-based limits have been presented recently |]l],|U]. 

In the next section we begin with a description of the data selection followed by a description of the calculation of 
backgrounds and observations of the data. Next we present rigorously-derived limits on both Minimal Supersymmetric 
(MSSM) and Gauge-mediated Supersymmetry Breaking (GMSB) Models. The next sections present the model- 
independent limits. Finally, in the Appendix we present tests of the application of model-independent limits to a 
variety of models that generate this signature. 

A search for the heavy Techniomega, ujt, in the final state 7 + 6 + jet, derived from the same data sample, has 
already been pubhshed p. 
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II. DATA SELECTION 



The data used here correspond to 85 pb^^ of pp cohisions at ^/s — 1.8 TeV. The data sample was collected by 
triggering on the electromagnetic cluster caused by the photon in the central calorimeter. We use 'standard' photon 
identification cuts developed for previous photon analyses which are similar to standard electron requirements 
except that there is a restriction on any tracks near the cluster. The events are required to have at least one jet with 
a secondary vertex found by the standard silicon detector 6-quark identification algorithm. Finally, we apply missing 
transverse energy requirements and other selections to examine subsamples. We discuss the selection in detail below. 

A. The CDF Detector 

We briefly describe the relevant aspects of the CDF detector ||^ . A superconducting solenoidal magnet provides 
a 1.4 T magnetic field in a volume 3 m in diameter and 5 m long, containing three tracking devices. Closest to 
the beamline is a 4-layer silicon microstrip detector (SVX) |Q] used to identify the secondary vertices from 6-hadron 
decays. A track reconstructed in the SVX has an impact parameter resolution of 19/im at high momentum to 
approximately 25/im at 2 GeV/c of track momentum. Outside the SVX, a time projection chamber (VTX) locates 
the z position of the interaction. In the region with radius from 30 cm to 132 cm, the central tracking chamber (CTC) 
measures charged-particle momenta. Surrounding the magnet coil is the electromagnetic calorimeter, which is in turn 
surrounded by the hadronic calorimeter. These calorimeters are constructed of towers, subtending 15° in (p and 0.1 in 
rj 1^, pointing to the interaction region. The central preradiator wire chamber (CPR) is located on the inner face of 
the calorimeter in the central region {\r]\ < 1.1). This device is used to determine if the origin of an electromagnetic 
shower from a photon was in the magnet coil. At a depth of six radiation lengths into the electromagnetic calorimeter 
(and 184 cm from the beamline), wire chambers with additional cathode strip readout (central electromagnetic strip 
chambers, CES) measure two orthogonal profiles of showers. 

For convenience we report all energies in GeV, all momenta as momentum times c in GeV, and all masses as mass 
times in GeV. Transverse energy (Et) is the energy deposited in the calorimeter multiplied by sin0. 

B. Event Selection 

Collisions that produce a photon candidate are selected by at least one of a pair of three-level triggers, each of which 
requires a central electromagnetic cluster. The dominant high-Et photon trigger requires a 23 GeV cluster with less 
than approximately 5 GeV additional energy in the region of the calorimeter surrounding the cluster |^ . A second 
trigger, designed to have high efficiency at large values of Et, requires a 50 GeV cluster, but has no requirement on 
the isolation energy. 

These events are required to have no energy deposited in the hadronic calorimeter outside of the time window that 
corresponds to the beam crossing. This rejects events where the electromagnetic cluster was caused by a cosmic ray 
muon which scatters and emits bremsstrahlung in the calorimeter. 

Primary vertices for the pp collisions are reconstructed in the VTX system. A primary vertex is selected as the 
one with the largest total \pt\ attached to it, followed by adding silicon tracks for greater precision. This vertex is 
required to be less than 60 cm from the center of the detector along the beamline, so that the jet is well-contained 
and the projective nature of the calorimeters is preserved. 

C. Photon 

To purify the photon sample in the offline analysis, we select events with an electromagnetic cluster with Et >25 GeV 
and 1 77 1 < 1.0. To provide for a reliable energy measurement we require the cluster to be away from cracks in the 
calorimeter. To remove backgrounds from jets and electrons, we require the electromagnetic cluster to be isolated. 
Specifically, we require that the shower shape in the CES chambers at shower maximum be consistent with that of 
a single photon, that there are no other clusters nearby in the CES, and that there is little energy in the hadronic 
calorimeter towers associated with (i.e. directly behind) the electromagnetic towers of the cluster. 

We allow no tracks with pt > 1 GeV to point at the cluster, and at most one track with pt < 1 GeV. We require 
that the sum of the pt of all tracks within a cone of AR = \/ Arf- -I- Ac^P- = 0.4 around the cluster be less than 5 GeV. 

If the photon cluster has Et < 50 GeV, we require the energy in a 3 x 3 array of trigger towers (trigger towers are 
made of two consecutive physical towers in rj) to be less than 4 GeV. This isolation energy sum excludes the energy 
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in the electromagnetic calorimeter trigger tower with the largest energy. This requirement is more restrictive than 
the hardware trigger isolation requirement, which is approximately 5 GeV on the same quantity. In some cases the 
photon shower leaks into adjacent towers and the leaked photon shower energy is included in the isolation energy 
sum. This effect leads to an approximately 20% inefhciency for this trigger. When the cluster Et is above 50 GeV, 
a second trigger with no isolation requirement accepts the event. For these events, we require the transverse energy 
found in the calorimeter in a cone of i? = 0.4 around the cluster to be less than 10% of the cluster's energy. 
These requirements yield a data sample of 511,335 events in an exposure of 85 pb^^ of integrated luminosity. 

D. B-quark Identification 

Jets in the events are clustered with a cone of 0.4 in 77 — space using the standard CDF algorithm ||l^ . One of the 
jets with Id < 2 is required to be identified as a 5-quark jet by the displaced-vertex algorithm used in the top-quark 
analysis |lj]. This algorithm searches for tracks in the SVX that are associated with the jet but not associated 
with the primary vertex, indicating they come from the decay of a long-lived particle. We require that the track, 
extrapolated to the interaction vertex, has a distance of closest approach greater than 2.5 times its uncertainty and 
pass loose requirements on pt and hit quality. The tracks passing these cuts are used to search for a vertex with three 
or more tracks. If no vertex is found, additional requirements are placed on the tracks, and this new list is used to 
search for a two-track vertex. The transverse decay length, L^y, is defined in the transverse plane as the projection 
of the vector pointing from the primary vertex to the secondary vertex on a unit vector along the jet axis. We require 
\Lxy\/<y > 3, where a is the uncertainty on L^y These requirements constitute a "tag". In the data sample the tag is 
required to be positive, with L^y > 0. The photon cluster can have tracks accidentally associated with it and could 
possibly be tagged; we remove these events. This selection reduces the dataset to 1487 events. 

The jet energies are corrected for calorimeter gaps and non-linear response, energy not contained in the jet cone, 
and underlying event energy [|o|. For each jet the resulting corrected Et is the best estimate of the underlying true 
quark or gluon transverse energy, and is used for all jet requirements in this analysis. We require the Et of the tagged 
jet in the initial 67 event selection to be greater than 30 GeV; this reduces the data set to 1175 events. 

E. Other Event Selection 

While the photon and 5-tagged jet constitute the core of the signature we investigate, supersymmetry and other 
new physics could be manifested in any number of different signatures. Because of the strong dependence of signature 
on the many parameters in supersymmetry, one signature is (arguably) not obviously more likely than any other. For 
these reasons we search for events with unusual properties such as very large missing Et or additional reconstructed 
objects. These objects may be jets, leptons, additional photons or 5-tags. This method of sifting events was employed 
in a previous analysis 0. We restrict ourselves to objects with large Et since this process is serving as a sieve of the 
events for obvious anomalies. In addition, in the lower Et regime the backgrounds are larger and more difficult to 
calculate. In this section we summarize the requirements that define these objects. 

Missing Et {f^t) is the magnitude of negative of the two-dimensional vector sum of the measured Et in each 
calorimeter tower with energy above a low threshold in the region \rj\ < 3.6. All jets in the event with uncorrected 
Et greater than 5 GeV and |r/| <2 are corrected appropriately for known systematic detector mismeasurements; these 
corrections are propagated into the missing Et. Missing Et is also corrected using the measured momentum of muons, 
which do not deposit much of their energy in the calorimeter. 

We apply a requirement of 20 GeV on missing Et, and observe that a common topology of the events is a photon 
opposite in azimuth from the missing Et (see Figure ||) . We conclude that a common source of missing Et occurs when 
the basic event topology is a photon recoiling against a jet. This topology is likely to be selected by the f^t cut because 
the fluctuations in the measurement of jet energy favor small jet energy over large. To remove this background, we 
remove events in the angular bin A<^(7— f^t) > 168° for the sample, where we have raised the missing Et requirement 
to 40 GeV. 

We define Ht as the scalar sum of the Et in the calorimeter added to the missing Et and the pt of any muons in 
the event. This would serve as a measure of the mass scale of new objects that might be produced. 

To be recognized as an additional jet in the event, a calorimeter cluster must have corrected Et >15 GeV and 
|?7| < 2. To count as an additional b tag, a jet must be identified as a 6 candidate by the same algorithm as the 
primary h jet, and have Et > 30 GeV and jr/j <2. To be counted as an additional photon, an electromagnetic cluster 
is required to have Et > 25 GeV, \rj\ < 1.0, and to pass all the same identification requirements as the primary photon. 
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Object 


Selection 


Basic Sample Requirements 


Isolated Photon 




Et > 25 GeV, \ri\ < 1.0 


fe-quark jet (SVX 6-tag) 




Et > 30 GeV, \ri\ < 2.0 


Optional Missing Et Requirements 






> 40 GeV 


|AV(7-^t)l 




< 168° 


Optional Other Objects 


Jets 




Et > 15 GeV, |r?| < 2.0 


Additional Photons 




Et > 25 GeV, \ri\ < 1.0 


Additional 6-quark jets 




Et > 30 GeV, SVX fe-tag 


Electrons 




Et > 25 GeV, \ri\ < 1.0 


Muons 




Pt > 25 GeV, \ri\ < 1.0 


Tau Leptons 




Et > 25 GeV, |r?| < 1.2 



TABLE I. Summary of the kinematic selection criteria for the + X sample that contains 1175 events. Also shown are 
the kinematic criteria for the identification of other objects, such as missing Et, jets, additional b-jets and leptons. The lepton 
identification criteria are the same as used in the top discovery. 



For lepton identification, we use the cuts defined for the primary leptons in the top quark searches |T^,[T^. We 
search for electrons in the central calorimeter and for muons in the central muon detectors. Candidates for r leptons 
are identified only by their hadronic decays - as a jet with one or three high-pt charged tracks, isolated from other 
tracks and with calorimeter energy cluster shapes consistent with the r hypothesis . Electrons and r's must have 
Et > 25 GeV as measured in the calorimeter; muons must havept > 25 GeV. Electrons and muons must have |7/| < 1.0 
while r's must have \r]\ < 1.2. We summarize the kinematic selections in Table |. 

III. BACKGROUND ESTIMATES 

The backgrounds to the Iry sample are combinations of standard model production of photons and b quarks and also 
jets misidentified as a photon ("fake" photons) or as a 6-quark jet ("fake" tags or mistags). A jet may be misidentified 
as a photon by fragmenting to a hard leading tt^. Other jets may fake a 6-quark jet through simple mismeasurement 
of the tracks leading to a false secondary vertex. 

We list these backgrounds and a few other smaller backgrounds in Table The methods referred to in this table 
are explained in the following sections. 



Source 


Method of Calculation 


"fbb and 7CC 


Monte Carlo 


7+ mistag 


CES-CPR and tagging prediction 


fake 7 and bb or cc 


CES-CPR 


fake 7 and a mistag 


CES-CPR 


W-y, Z-y 


Monte Carlo, normalized to data 


electrons faking 7's 


measured fake rate 


cosmic rays 


cosmic characteristics 



TABLE IL The summary of the backgrounds to the photon and tag sample and the methods used to calculate them. 
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The following sections begin with a discussion of the tools used to calculate backgrounds. Section IIIC explains 
why the method presented is necessary. The subsequent sections provide details of the calculation of each background 
in turn. 



A. Photon Background Tools 

There are two methods we use to calculate photon backgrounds, each used in a different energy region. The first 
employs the CES detector embedded at shower maximum in the central electromagnetic calorimeter . This method 
is based on the fact that the two adjacent photons from a high-pt tt" will tend to create a wide CES cluster, with 
a larger CES x^; when compared to the single photon expectation. The method produces an event-by-event weight 
based on the of the cluster and the respective probabilities to find this x^ for £^ ""^ versus for a photon. In the 
decay of very high-energy tt^'s the two photons will overlap, and the tt" will become indistinguishable from a single 
photon in the CES by the shape of the cluster. From studies of tt^'s from p decay we have found that for Et > 35 GeV 
the two photons coalesce and we must use a second method of discrimination that relies on the central preradiator 
system (CPR) [Q. This background estimator is based on the fact that the two photons from a 7r° have two chances 
to convert to an electron-positron pair at a radius before the CPR system, versus only one chance for a single photon. 
The charged particles from the conversion leave energy in the CPR, while an unconverted photon does not. The 
implementation of the CPR method of discriminating photons from tt'^'s on a statistical basis is similar to the CES 
method, an event-by-event weight. When the two methods are used together to cover the entire photon Et range for 
a sample, we refer to it as the CES-CPR method. 

Both these photon background methods have low discrimination power at high photon Et. This occurs because 
the weights for a single photon and a (background) 7r° are not very different. For example, in the CES method, at 
an Et of 25 GeV, the probability for a photon to have a large x^ is on the order of 20% while the background has 
a probability of approximately 45%. For the CPR method, typical values for a 25 GeV photon are 83% conversion 
probability for background and 60% for a single photon. 

B. 6-quark Tagging Background Tools 

A control sample of QCD multi-jet events is used to study the backgrounds to the identification of 6-quark jets [Q. 
For each jet in this sample, the Et of the jet, the number of SVX tracks associated with the jet, and the scalar sum 
of the Et in the event are recorded. The probability of tagging the jet is determined as a function of these variables 
for both positive {L^y > 0) and negative tags {L^y < 0). 

Negative tags occur due to measurement resolution and errors in reconstruction. Since these effects produce negative 
and positive tags with equal probability, the negative tagging probability can be used as the probability of finding a 
positive tag due to mismeasurement (mistags). 

C. Background Method 

We construct a total background estimate from summing the individual sources of backgrounds, each found by 
different methods. In the CDF top analysis one of the tagging background procedures was to apply the positive 
tagging probability to the jets in the untagged sample to arrive at a total tagging background estimate. A similar 
procedure could be considered for our sample. 

However, in this analysis, a more complex background calculation is necessary for two reasons. First, the param- 
eterized tagging background described above is derived from a sample of jets from QCD events pl| which have a 
different fraction of &-quark jets than do jets in a photon-plus-jets sample. This effect is caused by the coupling of the 
photon to the quark charge. Secondly, b quarks produce B mesons which have a large branching ratio to semileptonic 
states that include neutrinos, producing real missing Et more often than generic jets. When a ^t cut is applied, the b 
fraction tends to increase. This effect is averaged over in the positive background parameterization so the background 
prediction will tend to be high at small and low at large 

For these reasons, the positive tagging rate is correlated to the existence of a photon and also the missing Et, when 
that is required. In contrast, the negative tagging rate is found not to be significantly correlated with the presence 
of real b quarks. This is because the negative tagging rate is due only to mismeasurement of charged tracks which 
should not be sensitive to the flavor of the quarks. 
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The next sections list the details of the calculations of the individual sources of the backgrounds. Both photons and 
6-tagged jets have significant backgrounds so we consider sources with real photons and &-tags or jets misidentified 
as photons or &-jets ("fakes"). 



D. Heavy Flavor Monte Carlo 



The background consisting of correctly-identified photons and 6-quark jets is computed with an absolutely normal- 
ized Monte Carlo The calculation is leading order, based on qq and gg initial states and a finite 6-quark mass. 
The scale is taken to be the square of the photon Et plus the square of the bb or cc pair mass, = E^^ + M^. A 
systematic uncertainty of 30% is found by scaling Q by a factor of two and the quark masses by 10%. An additional 
20% uncertainty allows for additional effects which cannot be determined by simply changing the scale dependence 
@- 

In addition, we rely on the detector simulation of the Monte Carlo to predict the tail of the rapidly falling 
spectrum. The Monte Carlo does not always predict this tail well. For example, a Monte Carlo of Z ^ e+e~ 
production predicts only half the observed rate for events passing the missing Et cut used in this analysis. We thus 
include an uncertainty of 100% on the rate that events in the &7 sample pass the Pt cut. We combine the Monte Carlo 
production and sources of uncertainty in quadrature. However when the 766 and 7CC backgrounds are totaled, 
these common uncertainties are treated as completely correlated. 



E. Fake photons 



The total of all backgrounds with fake photons can be measured using the CES and CPR detectors. These 
backgrounds, dominated by jets that fragment to an energetic tt" 77 and consequently are misidentified as a 
single photon, are measured using the shower shape in the CES system for photon Et < 35 GeV and the probability 
of a conversion before the CPR for Et > 35 GeV [|l6|. We find 55 ± 1 ± 15% ||l^ of these photon candidates are 
actually jets misidentified as photons. 

For many of our subsamples we find this method is not useful due to the large statistical dilution as explained in 



Section This occurs because, for example, the probabilities for background (tt^'s) and for signal (7's) to convert 
before the CPR are not too different. This results in a weak separation and a poor statistical uncertainty. We find 
the method returns 100% statistical uncertainties for samples of less than approximately 25 photon candidates. 



F. Real photon, Fake tags 



To estimate this background we start with the untagged sample, and weight it with both the CES-CPR real photon 
weight and the negative tagging (background) weight. This results in the number of true photons with mistags in the 
final sample. As discussed above, the negative tagging prediction does not have the correlation to quark flavor and 
missing Et as does the positive tagging prediction. 

As a check, we can look at the sample before the tagging and requirements. In this sample we find 197 negative 
tags while the estimate from the negative tagging prediction is 312. This discrepancy could be due to the topology of 
the events - unlike generic jets, the photon provides no tracks to help define the primary vertex. The primary vertex 
could be systematically mismeasured leading to mismeasurement of the transverse decay length L^y for some events. 
We include a 50% uncertainty on this background due to this effect. 



G. Estimate of Remaining Backgrounds 



There are several additional backgrounds which we have calculated and found to very small. The production of 
and Zj events may provide background events since they produce real photons and 6 or c quarks from the boson 
decay (W^ ^ cs, Z ^ bb). The f]^ would have to be fake, due to mismeasurement in the calorimeter. We find W/Zj 
events in the CDF data using the same photon requirements as the search. The W/Z is required to decay leptonically 
for good identification. We then use a Monte Carlo to measure the ratio of the number of these events to the number 
of events passing the full 76^4 search cuts. The product of these two numbers predicts this background to be less 
than 0.1 events. 
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The next small background is W ^ eiy plus jets where the electron track is not reconstructed, due either to 
bremsstrahlung or to pattern-recognition failure. Using Z — > e'^e~ events, we find this probability is small, about 
0.5%. Applying this rate to the number of observed events with an electron, fe-tag and missing Et we find the number 
of events expected in our sample to be negligible. 

The last small background calculation is the rate for cosmic ray events. In this case there would have to be a QCD 
&-quark event with a cosmic ray bremsstrahlung in time with the event. The missing Et comes with the unbalanced 
energy deposited by the cosmic ray. We use the probability that a cosmic ray leaves an unattached stub in the muon 
chambers to estimate that the number of events in this category is also negligible. 



The total of all background sources is summarized in Table III , The number of observed events is consistent with 



the calculation of the background for both the 76 sample and the subsamples with 



IV. DATA OBSERVATIONS 



In this section we report the results of applying the final event selection to the data. First we compare the 
total background estimate with the observed number of events in the 67 sample, which requires only a photon with 
Et > 25 GeV and a 6-tagged jet with Et > 30 GeV. Since most models of supersymmetry predict missing Et, we also 
tabulat e th e backgrounds for that subsample. 

Table HH summarizes the data samples and the predicted backgrounds. We find 98 events have missing Et > 20 GeV. 
Six events have missing Et > 40 GeV, and only two of those events pass the A(/)(7 — ^f ) < 168° cut. 

Figures and || display the kinematics of the data with a background prediction overlayed. Because of the large 
statistical uncertainty in the fake photon background, the prediction for bins with small statistics have such large 
uncertainties that they are not useful. In this case we approximate the fake photon background by applying the fake 
photon measurement and the positive tagging prediction to the large-statistics untagged sample. This approximation 
only assumes that real b quarks do not produce substantial missing Et. Each component of the background is 



normalized to the number expected as shown in Table 111; the total is then normalized to the data in order to 



compare distributions. We observe no significant deviations from the expected background. 



Source 


Events 


Events > 20 


Events fJ^ > 40, 


766 


99 ± 5 ± 50 


9± 1 ± 10 


0.4 ±0.3 ±0.4 


7CC 


161±9±81 


7±2±8 


0.0 ±0.5 ±0.5 


7-l-mistag 


124 ± 1 ± 62 


10 ±0.3 ±5.2 


0.7 ±0.05 ±0.5 


fake 7 


648 ± 69 ± 94 


49 ± 22 ± 7 


1.0 ±1.0 ±0.2 


W7 


2± 1 


0.4 ±0.2 ±0.4 


0.0 ±0.1 ±0.1 


Zl 


6±4 


0.8 ±0.6 ±0.8 


0.08 ±0.06 ±0.08 


e ^ 7 


0.4 ±0.1 


0.4 ±0.1 


0.1 ± .03 


cosmics 


0± 16 


0±5 





total background 


1040 ± 72 ± 172 


77 ± 23 ± 20 


2.3 ±1.2 ±1.1 


data 


1175 


98 


2 



TABLE III. Summary of the primary background calculation. The 766 and 7CC systematic uncertainties are considered 
100% correlated. The column labeled > 40 GeV also includes the requirement that A(j){'y — p^) < 168°. The entry for fake 
photons in the column labeled > 40 GeV is not measured but is estimated using the assumption that 50% of photons are 
fakes. This number is assigned a 100% uncertainty. 



Several events appear on the tails of some of the distributions. Since new physics, when it first appears, will likely 
be at the limit of our kinematic sensitivity, the tail of any kinematic distribution is a reasonable place to look for 
anomalous events. However, a few events at the kinematic limit do not warrant much interest unless they have many 
characteristics in common or they have additional unusual properties. We find two events pass the largest missing 
Et cut of 40 GeV; we examine those events more closely below. We also obser ve th ere are five events with large dijet 
mass combinations and we also look at those more closely below. In Section [V C , we search for other anomalies in 
our sample. 
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FIG. 1. Comparison of the data to the background prediction (dashed hne), and the basehne SUSY model of Section V A2 
(dotted line). The data consist of the 98 events of the 76 data with > 20 GeV, except in b) which contains no 
requirement. In each case the predictions have been normalized to the data. The distributions are: a) the photon Et, b) the 
missing Et, c) the fe-tagged jet Et and d) the Et of the second jet with Et >15 GeV, if there is one. For display, the SUSY 
model event yield is scaled up by a factor of 4 for a), c) and d) and a factor of 40 for b). 
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FIG . 2. Comparison of the data to the background prediction (dashed line), and the the baseline SUSY model of Section 
V A ^ (dotted line), each normalized to the 98 events of the 76 data with fi^ > 20 GeV. The distributions are: a) Ht (total 
energy), b) between the photon and the p^, c) number of jets with Et > 15 GeV, and d) A(j) between the missing Et and 
the nearest jet. For display, the SUSY model event yield is scaled up by a factor of 4. 
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FIG. 4. The distributions for: a) M{b,j) and b) M{'y,b,j) for the fi^ > 20 GeV events as shown in Figure ||. Only 63 of the 
98 events have a second jet and ma ke it in to this plot. The data are compared to a background prediction (dashed line), and 
the baseline SUSY model of Section V A2 (dotted line) , each normalized to the data. The Monte Carlo prediction is scaled up 
by a factor of 3. 
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A. Analysis of Events with Large Missing Et 



Six events pass the a priori selection criteria requiring a photon, 6-tag, and >40 GeV. Two of these events also 
pass the A(/)(7 — < 168° requirement. We have examined these two events to see if there indications of anything 
else unusual about them (for example, a high-pt lepton, or a second jet which forms a large invariant mass with the 
first 5-jet, to take signals of GMSB and Higgsino models respectively). 

The first event (67537/59517) does not have the characteristics of a typical 5-tag. It is a two-track tag (which has 
a worse signal-to-noise) with the secondary vertex consistent with the beam pipe radius (typical of an interaction 
in the beam pipe). The two tracks have a pt of 2 and 60 GeV, respectively; this highly asymmetric configuration is 
unlikely if the source is a 6-jet. There are several other tracks at the same 4> as the jet that are inconsistent with 
either the primary or secondary vertex. We conclude the 6-tag jet in this event is most likely to be a fake, coming 
from an interaction in the beam pipe. 

The second event has a typical 6-tag but there are three jets, and all three straddle cracks in the calorimeter 
(77 = 0.97, —1.19, —0.09), implying the is very likely to be mismeasured. 

In both events we judge by scanning that the primary vertex is the correct choice so that a mismeasurement of the 
^( due to selecting the wrong vertex is unlikely. 

While we have scanned these two events and find they are most likely not true jb^^ events, we do not exclude 
them from the event sample as the background calculations include these sources of mismeasured events. 



Run/Event 


7 Et 




M(b,jet) 


b Et 


jets Et 


A<^(7-^t) 


A(j}{j„ear ~ ^t) 


Ht 


60951/189718 


121 


42 


57 


61 


67,26,15 


177 


11 


342 


64997/119085 


222 


44 


97 


173 


47 


170 


1 


495 


63684/15166 


140 


57 


63 


35 


25,20,15 


175 


6 


388 


67537/59517* 


36 


73 


399 


195 


141,113,46,17 


124 


20 


595 


69426/104696 


33 


58 


266 


143 


119 


180 


3 


344 


68464/291827* 


93 


57 


467 


128 


155,69 


139 


16 


405 



TABLE IV. Characteristics of the six events with > 40 GeV; the two marked with an asterisk also pass the 
A(f){'y — < 168° requirement. All units are GeV except for which is in degrees. The columns are the Et of the 
photon in the event, the missing Et, the mass of the 6-jet and the second highest Et jet in the event, the Et of the b jet, the Et 
of the jets other than the b jet, the Acp between the photon and the missing Et, the A<j) between the the missing Et and the 
nearest jet, and the Ht of the event (scalar sum of the Et, the missing Et and the pt of any muons in the event). 



B. Analysis of Five High— mass Events 

If the events include production of new, heavy particles, we might observe peaks, or more likely, distortions in the 
distributions of the masses formed from combinations of objects. To investigate this, we create a scatter plot of the 
mass of the 6-quark jet and the second highest-£'t jet versus the mass of the photon, &-quark jet and second highest 
jet in Figure || and |^. 

As seen in the figures, the five events at highest M{b, j) seem to form a cluster on the tail of the distribution. There 
are 63 events in the scatter plot which are the subset of the 98 events with > 20 GeV which contain a second jet. 
The five events include the two (probable background) events with > 40 GeV and A0(7 — p^) < 168° and three 
events with large Ht (> 400 GeV). Since these events were selected for their high mass, we expect they would appear 
in the tails of several of the distributions such as Ht. Table M shows the characteristics of these five events. 



Run/Event 


7 Et 


f^t 


M{b,jet) 


b Et 


jet Ets 


A<^(7 - ^t) 


A(j){j„ear - ^ t) 


Ht 


66103/52684 


106 


24 


433 


170 


135,57 


152 


29 


517 


66347/373704 


122 


32 


369 


268 


125,42 


101 


14 


605 


67537/59517 


36 


73 


399 


195 


141,113,46,17 


124 


20 


595 


68333/233128 


38 


39 


395 


99 


282,212 


121 


3 


600 


68464/291827 


93 


57 


467 


128 


155,69 


139 


16 


405 



TABLE V. Characteristics of events with M{b,jet) > 300 GeV. For a complete description of the quantities, see Table 
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In order to see if these events are significant, we need to make an estimate of the expected background. We define 
the two regions indicated in Figure ||. The smaU box is placed so that it is close to the five events]^ This is intended 
to maximize the significance of the excess. The large box is placed so that it is as far from the five events as possible 
without including any more data events. This will minimize the significance. The two boxes can serve as informal 
upper and lower bounds on the significance. Since these regions were chosen based on the data, the excess over 
background cannot be used to prove the significance of these events. These estimates are intended only to give a 
guideline for the significance. 

We cannot estimate the background to these five events using the CES and CPR methods described in Section [II A 
due to the large inherent statistical uncertainties in these techniques. We instead use the following procedure. The 
list of backgrounds in Section [II defines the number of events from each source with no restriction on M{b,j). We 
normalize these numbers to the 63 events in the scatter plot. We next derive the fraction of each of these sources we 
expect at high AI{b,j). We multiply the background estimates by the fractions. The result is a background estimate 
for the high-mass regions. 

To derive the fractions of background sources expected at high M{b,j) we look at each background in turn. The 
fake photons are QCD events where a jet has fluctuated into mostly electromagnetic energy. For this source we use 
the positive Lxy background prediction [ pT| to provide the fraction. This prediction is derived from a QCD jet sample 
by parameterizing the positive tagging probability as a function of several jet variables. The probability for each jet 
is summed over all jets for the untagged sample to arrive at a tagging prediction. Since the prediction is derived from 
QCD jets we expect it to be reliable for these QCD jets also. Running this algorithm (called "Method 1" ||ll|]) on the 
untagged photon and ff^ sample yields the fraction of expected events in each of the two boxes. The fractions are 



summarized in Table Vl 



Source 


Big box 


Small box 


fake 7 


0.080 ± 0.007 


0.017 ±0.003 


7, fake tag 


0.112 ±0.009 


0.032 ± 0.005 


fbb 


0.10 ±0.03 


0.022 ± 0.007 


^cc 


0.08 ± 0.04 


0.018 ± 0.008 



TABLE VI. The fraction of the 63 ybj f]^ events for each background expected to fall into the high-A/(6, j) boxes defined in 
Figure ^ 



^Note that events cannot be above the diagonal in the M(7, b,j) — M{b,j) plane, so the true physical area is triangular. 
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The second background source considered consists of real photons with fake tags. We calculate this contribution 
using the measured negative tagging rate applied to all jets (i.e. before 6-tagging) in the sample. Finally, the real 
photon and heavy flavor backgrounds are calculated based on the Monte Carlo. The results from estimating the 
fractions are shown in Table 



VI 



The estimates of the sources of background for the 63 events at all M{b,j) have statistical uncertainties, as do 



the fractions in Table VI; we include both in the uncertainty in the number of events in the high-mass boxes. We 
propagate the systematic uncertainties on the backgrounds to the 63 events at all M{b,j) and include the following 
systematics due to the fractions: 

1. 50% of the real photon and mistag background calculation for the possibility that the quark and giuon content, 
as well as the heavy flavor fraction, in photon events may differ from the content in QCD jets. 

2. 50% of the real photon and mistag background calculation for the possibility that using the positive tagging 
prediction to correct the Monte Carlo for the cut may have a bias. 

3. 100% of the real photon and real heavy flavor background calculation for the possibility that the tails in the 
Monte Carlo M (6, j) distribution may not be reliable. 

The results of multiplying the backgrounds at all M{b,j) with the fractions expected at high M{b,j) are shown in 



Table VII. 



Source 


Big box 


Small box 


fake 7 


3.3 ±1.5 ±0.5 


0.70 ±0.33 ±0.10 


7, fake tag 


0.97 ±0.09 ±0.69 


0.28 ±0.05 ±0.20 


ybb 


0.75 ±0.26 ±1.18 


0.16 ±0.06 ±0.26 


ICC 


0.44 ±0.26 ±0.79 


0.11 ±0.06 ±0.17 


total 


5.5 ±1.5 ±1.6 


1.24 ±0.35 ±0.38 



TABLE VII. Summary of the estimates of the background at high M{b,j) in the boxes in the M(7,6,j) — M{b,j) plane 
defined in Figure B 



The result is that we expect 5.5 ± 1.5 ± 1.6 events in the big box, completely consistent with five observed. We 
expect 1.2 ± 0.35 ± 0.38 events in the small box. The probability of observing five in the small box is 1.6%, a 2.7(7 
effect, a posteriori. 

We next address a method for avoiding the bias in deciding where to place a cut when estimating backgrounds 
to events on the tail of a distribution. This method was developed by the Zeus collaboration for the analysis of the 
significance of the tail of the distribution [Q. Figure ^ summarizes this method. The Poisson probability that 
the background fluctuated to the observed number of events (including uncertainties on the background estimate) is 
plotted for different cut values. We use the projection of the scatter plot onto the M{b,j) axis and make the cut on 
this variable since this is where the effect is largest. We find the minimum probability is 1.4 x 10""^, which occurs 
for a cut of M{b,j) >400 GeV. We then perform 10,000 "pseudo-experiments" where we draw the data according 
to the background distribution derived above and find the minimum probability each time. We find 1.2% of these 
experiments have a minimum probability lower than the data, corresponding to a 2.7a fiuctuation. Including the 
effect of the uncertainties in the the background estimate does not significantly change the answer. 

We note that this method is one way of avoiding the bias from deciding in what region to compare data and 
backgrounds after seeing the data distributions. It does not, however, remove the bias from the fact that we are 
investigating this plot, over all others, because it looks potentially inconsistent with the background. If we make 
enough plots one of them will have a noticeable fiuctuation. We conclude that the five events on the tail represent 
something less than a 2.7cr effect. 
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FIG. 5. M{b,j) versus M{'y,b,j) for the events with > 20 GeV as shown in Figure ^ Only 63 of the 98 events have a 
second jet and make it into this plot. The small dots are the result of making the scatter plot for the untagged data (passing 
all other cuts) and weighting it with the positive tagging prediction. The estimates of background expected in the boxes are 
found by the method described in the text. 
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FIG. 6. The upper plot is the number of events passing a cut on M{b, j) for the data and the positive tagging prediction. The 
lower plot is the probability that the number of events passing a cut on M{b, j) is consistent with the positive tagging prediction. 
The expected number of experiments with such a low minimum probability is derived from 10,000 simulated experiments drawn 
from the distribution of the expected background. 
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C. Additional Objects in the Data Sample 



We have searched the 76 data sample for other unusual characteristics. The creation and decay of heavy squarks, 
for example, could produce an excess of events with multiple jets. In Figure |^ we histogram the number of events 
with N or more jets. Table VIII presents the numbers of events observed and expected. Some backgrounds are 
negative due to the large statistical fluctuations of the fake photon background. When all backgrounds are included 
the distribution in the number of jets in the data is consistent with that from background. 



Min Njet 


Observed, > GeV 


Expected, > GeV 


Observed, > 20 GeV 


Expected, > 20 GeV 


1 


1175 


1040 ± 72 ± 172 


98 


77 ± 23 ± 20 


2 


464 


394 ± 44 ± 63 


63 


39 ±18 ±12 


3 


144 


82 ± 24 ± 14 


25 


-8± 12±3 


4 


36 


17± 11 ±3 


7 




5 


10 




3 




6 


5 




1 




7 


2 









8 


1 










TABLE VIII. Numbers of events with A'^ or more jets and the expected Standard Model background. Some background 
predictions are negative due to the large statistical fluctuations on the fake photon background method. 



We have searched in the sam ple o f events with a photon and ^-tagged jet for additional high-iJt objects using the 
requirements defined in Section [IE. We find no events containing a second photon. We find no events containing a 
hadronic r decay or a muon. We find one event with an electron; its characteristics are listed in Table IX. In scanning 
this event, we note nothing else unusual about it. 

We find 8 events of the 1175 which have a photon and 5-tagged jet contain a second 6-tagged jet with Et >30 GeV. 
(Out of the 1175, only 200 events have a second jet with Et >30 GeV.) Unfortunately, this is such a small sample 
that we cannot use the background calculation to find the expected number of these events (the photon background 
CES~CPR method returns 100% statistical uncertainties). One of the events with two tags has 30 GeV of missing Et 
so is in the 98~event > 20 GeV sample. 



Run/Event 


7 Et 


fit 


M(7,e) 


b Et 


electron Et 


A<^(7 - fit) 


Ht 


63149/4148 


42 


17 


21 


106 


33 


43 


212 



TABLE IX. Characteristics of the one event with a photon, tagged jet, and an electron. 
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FIG. 7. The distribution of the number of events with A*' or more jets, represented by the sohd points. The boxes are centered 
on the background prediction and their size reflects plus and minus one-cr of combined statistical and systematic uncertainty 
on the background prediction. The distributions are: a) all events with a photon and 6-tagged jet, b) all events with a photon, 
6-tagged jet and > 20 GeV. Some background predictions are negative due to the large statistical fluctuations on the fake 
photon background method. The results are also tabulated in Table VIII. 
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V. LIMITS ON MODELS OF SUPERSYMMETRY 



In the following sections we present limits on three specific models of supersymmetry [ p^ . Each of these models 
predicts significant numbers of events with a photon, a 5-quark jet and missing transverse energy (ie 76^4). 

As is typical for supersymmetry models, each of these shows the problems in the process of choosing a model and 
presenting limits on it. Each of these models is very specific and thus represents a very small area in a very large 
parameter space. Consequently the odds that any of these is the correct picture of nature is small. They are current 
theories devised to address current concerns and may appear dated in the future. (This aspect is particularly relevant 
to the experimentalists, who often publish their data simultaneously with an analysis depending on a current model.) 
In addition these models can show sensistivity to small changes in the parameters. 

The first model is based on a particular location in MSSM parameter space which produces the signature of 'yb^^+X. 
We consider both direct production of charginos and neutralinos and, as a second model, indirect production of 
charginos and neutralinos through squarks and gluinos. The third model is based on the gauge-mediated concept, 
discussed further below. 



A. X2 ^ 1X1 Model 



This theoretical model was originally proposed in the context of the anomalous CDF 6677 event y,^. Here, 
however, we go beyond the constraints of this single event and only retain the essential elements of the model, 
optimized for CDF detector acceptence and efficiency. This is an MSSM model without any specific relation to a 
high-energy theory. It does not assume high-energy constraints such as the unification of the sfermion or scalar 
masses as is assumed in the models inspired by Supergravity (SUGRA) jl9j . In this section and the following section 
we develop a baseline model point in parameter space. The final limits on this model will be found for this point and 
for some variations around this point. 



1. Direct Gaugino Production in X2 ~* 7Xi Model 

The first part of the model ||] is a light stop squark (f), the superpartner to the top quark. In this model the xf 
then decays to tb and the t decays to XiC The second important feature of the model is the decay X2 — > 7X1 which 
dominates in a particular region of MSSM parameter space. With these decays dominating, any event where a xf 
and a X2 is produced, either directly or indirectly through the strong production and decays of squarks and gluinos, 
will contain a photon, a &-quark jet and missing Et- 

The heart of the model g] is the decay X2 ~^ 7Xi so we examine in detail the parameter space where this decay 
dominates. The branching ratio of this decay is large when one of the neutralinos is a pure photino and one is pure 
Higgsino. To make a pure photino, we set Mi = M2. The photino mass is then equal to M2. To make a pure Higgsino 
we set tan/3 « 1. To avoid the theoretical bias against a very small tan/3 (which makes the top Yukawa coupling go 
to infinity before the GUT scale) we will use tan/3 = 1.2. In this case, the Higgsino mass is approximately equal to 
\fj,\. The above is purely a result of the form of the neutralino mass matrix. For definitions of these model parameters 
and discussions of their roles in SUSY models please see jlH . 

This leaves two free parameters to define the charginos and neutralinos, M2 and /i. Figure |^ shows five regions in 
the n — M2 plane; Table ^ summarizes the regions. First we note that in region 5 (/i > 0) we do not observe the 
decay X2 ~^ 7Xi because typically < X2 X2 ~^ W*Xi- 

For /i < 0, there are four regions. In region 2, which is the region suggested in ||^, the X2 is the photino, Xi is the 
Higgsino, and the decay X2 ~* 7Xi dominates. In region 3, X2 is the Higgsino and Xi is the photino and the photon 
decay still dominates. In region 1 the photino has become so heavy it is now the Xs- In region 4, the Higgsino has 
become the Xs- In regions 3 and 4 it is still possible to get photon decays, sometimes even X3 7X2- 



Region 


X? 




xi 


xl 


1 


hb 


z 


7 


ha 


2 


hb 


7 


Z 


ha 


3 


7 


hb 


z 


ha 


4 


7 


Z 


hb 


ha 



TABLE X. The approximate content of the neutralinos in the four regions of the /i — A/2 plane with 11 < shown in Figure 
^ The symbols ha and hb are simply the antisymmetric and symmetric combinations of Hi and H2 [H . 



22 



We choose to concentrate on region 2 where the photon plus b decay signature can be rehably estimated by the 
Monte Carlo event generator PYTHIA The X2 ~^ 7Xi decay dominates here. We also note that in this region the 
cross section for X2X2 is 3-10 times larger than the cross section for xfx2 ^ven though the xt significantly heavier 
than the xf ■ This is due to the large W component of the x^ . 

Since region 2 is approximately one-dimensional, we scan in only one dimension, along the diagonal, when setting 
limits on X2X2 production. To decide where in the region to place the model, we note that the mass of X2 equals 
M2 and the mass of Xi = ImI in this region. To give the photon added boost for a greater sensitivity, we will set M2 
significantly larger than This restricts us to the upper part of region 2. The dotted line in Figure ^ is the set of 
points defined by all these criteria and is given by M2 — 0.89 * + 39 GeV. 

The next step is to choose a t mass. It is necessary that Xi < t < xt ^'^^ the decay Xi^:X2^ ~^ bi to dominate. 
We find that in Region 2, xt ~ ^2- If the i mass is near the xf , the b will only have a small boost, but the Xi in 
the decay t — > ex? will have a greater boost, giving greater If the i mass is near the Xi) the opposite occurs. 
In Monte Carlo studies, we find considerably more sensitivity if the t mass is near the Xi- We set the t mass to be 
M^o +5 GeV. Since the production cross section is larger than xfx^ ^'^d will be detected with better efficiency, 

when we simulate direct production we set the Monte Carlo program to produce only X^xH pairs. The final limit 
is expressed as a cross section limit plotted versus the X2 mass (which is very similar to the X2 mass). This model 
is designed to provide a simple, intuitive signature that is not complicated by branching ratios and many modes of 
production. 

For the baseline model, we chose a value of n near the exclusion boundary of current limits ||2l[| on a t which decays 
to cXi- The point we chose is M^o = 80 GeV. From the above prescription, this corresponds to Af^o = — = 80 GeV, 
M^o = M^± = M2 = 110 GeV, and = 85 GeV. This point, indicated by the dot in Figure g, gives the hghtest 
mass spectrum with good mass splittings that is also near the exclusion boundary from LEP and DO Collaborations. 

2. Squarks and Gluinos 

Now we address the squarks and gluinos, which can produce xfx2 in their decays, and sleptons, which can appear 
in the decays of charginos and neutralinos. 

We will set the squarks (the lighter b and both left and right u, d, s and c) to 200 GeV and the gluino to 210 GeV. 
The heavier t and b are above 1 TeV. The gluino will decay to the squarks and their respective quarks. The squarks will 
decay to charginos or neutralinos and jets. This will maximize the production of xfx2 and therefore the sensitivity. 

This brings us to the limit on indirect production in the X2 ~^ 7Xi model. The chargino and neutralino parameters 
are fixed at the baseline model parameters. We then vary the gluino mass and set the squark mass according to 
Mg = Mq + 10 GeV. The limit is presented as a limit on cross section plotted versus the gluino mass. When 

the gluino mass crosses the it threshold at 260 GeV, the gluino can decay to it and production of xtx2 decreases. 
However, since all squarks are lighter than the gluino, the branching ratio to the t is limited and production will not 
fall dramatically. 

Some remaining parameters of the model are now addressed. Sleptons could play a role in this model. They have 
small cross sections so they are not often directly produced, but if the sleptons are lighter than the charginos, the 
charginos can decay into the sleptons. In particular, the chargino decay tb may be strongly suppressed if it competes 
with a slepton decay. We therefore set the sleptons to be very heavy so they do not compete for branching ratios. We 
set Ma large. The lightest Higgs turns out to be only 87 GeV due to the corrections from the light third-generation 
squarks. This is below current limits so we attempted to tune the mass to be heavier and found it was difhcult to 
achieve, given the light t and low tan /3. 

Using the PYTHIA Monte Carlo program, we find that 69% of all events generated with squarks and gluinos have 
the decay X2 ~^ 7Xi; 58% have the decay xf ib, and 30% have both. (To be precise, the light stop squark was 
excluded from this exercise, as it decays only to cxi- A light stop pair thus gives the signature cc + one of the 
signatures used to search for it, |^,^ but not of interest here.) 
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3. Acceptances and Efficiencies 



This section describes the evaluation of the acceptance and efficiency for the indirect production of xfXi through 
squarks and gluinos and the direct production of X2X2 ^^e MSSM model of X2 ~^ 7Xi- We use the PYTHIA Monte 
Carlo with the CTEQ4L parton distribution functions (PDFs) 1 23 1 . The efficiencies for squark and gluino production 
at the baseline point are listed in Table Xl. 



Cuts 




Cumulative Efficiency (%) 


Photon 


Et > 25 GeV,|r7| < 1.0, ID cuts 


50 


One jet 


Et.corr > 30 GeV, \v\ < 2.0 


47 


One SVX tag 


Et.corr > 30 GeV, \ri\ < 2.0 


4.3 




> 40 GeV 


2.9 



TABLE XI. Efficiencies for the baseline point with squark and gluino production. The efficiencies do not include branching 
ratios. 
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The total efficiencies, which will be used to set p rodu ction limits below, are listed in Table XII for the production 
F xfxi through squarks and glui] 
former case, and 1% in the latter. 



of xfXi through squarks and gluinos, and in Table Kill for direction production. Typical efficiencies are 2-3% in the 



4- Systematic Uncertainty 

Some systematics are common to the indirect production and the direct production. The efficiencies of the isolation 
requirement in the Monte Carlo and Z — s- e^e~ control sample cannot be compared directly due to differences in the 
i?t-spectra of the electromagnetic cluster, and the multiplicity and Et spectra of associated jets. The difference (14%) 
is taken to be the uncertainty in the efficiency of the photon identification cuts. The systematic uncertainty on the 
^-tagging efficiency (9%) is the statistical uncertainty in comparisons of the Monte Carlo and data. The systematic 
uncertainty on the luminosity (4%) reflects the stability of luminosity measurements. 

We next evaluate systematics specifically for the indirect production. The baseline parton distribution function is 
CTEQ4L. Comparing the efficiency with this PDF to the efficiencies obtained with MRSDO' [|| and GRV~94LO H] 
for the squark and gluino production, we find a standard deviation of 5%. Turning off initial- and final-state radiation 
(ISR/FSR) in the Monte Carlo increases the efficiency by 1% and 2% respectively and we take half of these as the 
respective systematics. Varying the jet energy scale by 10% causes the efficiency to change by 4%. In quadrature, the 
total systematic for the indirect production is 18%. 

Evaluating the same systematics for the direct production, we find the uncertainty from the choice of PDF is 5%, 
from ISR/FSR is 2%/9%, and from jet energy scale is 4%. In quadrature, the total systematic uncertainty for the 
direct production is 20%. 



5. Limits on the xl ^ 7Xi Model, Indirect Production 



To calculate an approximate upper limit on the number oi^bp^ events from squark and gluino production, we use 
the limit implied from the observed 2 events, including the effect of the systematic uncertainties [^|7|. We divide 
the Poisson probability for observing <2 events for a given expected signal and background, convoluted with the 
uncertainties, by the Poisson probability for observing <2 events for a given background only, also convoluted with 
the uncertainties. The number of expected signal events is increased until the ratio falls below 5%, leading to an 
approximate 95% confidence level limit of 6.3 events. Other limits in this paper are computed similarly. 

This upper limit, the efficiency described above (also see Table XII), and the luminosity, 85 pb^^, are combined 
to find the cross section limit for this model. The theoretical cross section is calculated at NLO using the PRDSPIND 
program [E8[. The effect is to uniformly increase the strong interaction production cross sections by 30% (improving 
the limit). T_t the baseline point (including squarks and gluinos) described above, we expect 18.5 events, so this point 
is excluded. Next we find the limit as a function of the gluino mass. The squark mass is 10 GeV below the gluino 
mass and the rest of the sparticles are as in t he b aseline point. We can exclude gluinos out to a mass of 245 GeV in 
this model. The limits are displayed in Table XII and Figure 0. 



MgiGeV) 


Mij (GeV) , 


Gth X BR (pb) 


Ae (%) 


0-95% Hm X BR (pb) 


185 


175 


16.8 


1.97 


3.76 


210 


200 


7.25 


2.98 


2.49 


235 


225 


3.49 


3.23 


2.30 


260 


250 


1.94 


2.69 


2.76 


285 


275 


1.24 


2.16 


3.45 



TABLE XII. Efficiency times acceptance and limits on indirect production of X^X2 though squarks and gluinos. Approx- 
imately 30% of events contain the decays X2 ~* 7Xi smd ib. The efficiencies in this table are found as the number 
of events passing all cuts divided by the number of events that contain both of these decays. The product of cross section 
times branching ratio listed in each case is for all open channels of SUSY production. Masses are given in GeV (following our 
convention of quoting M x c^) and cross sections are in pb. The second row is the baseline point. 
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6. Limits on the X2 — > 7Xi Model, Direct Production 



In this case the number of observed events (two) is convoluted with the systematic uncertainty to obtain an upper 
Hmit of 6.4 events (95% C.L.). To calculate the expected number of events from the direct production of x^X2j we 



10, For 



vary /i, and calculate the and M2 as prescribed above. The results are shown in Table XII] and Figure 
these values of the model parameters, the branching ratios lb and X2 ~^ 7Xi s-re 100%^ 

As can be seen from Figure [Ifl, the predicted rates from direct production are smaller than the measured limits by 



1-3 orders-of-magnitude, and no mass limit on the X2 ma-ss can be set. 



M-o = -M 








Mi 


O'th 


Ae (%) 


0"95% Urn 


25 


61 


71 


110 


30 


0.23 


0.93 


8.06 


62 


95 


94 


130 


67 


0.034 


1.41 


5.33 


79 


110 


108 


140 


85 


0.018 


1.29 


5.85 


93 


123 


118 


150 


98 


0.0075 


1.34 


5.58 


118 


146 


140 


170 


123 


0.0022 


1.27 


5.94 



TABLE XIII. Efficiencies and limits on direct production of xjx2- Branching ratios xi ~^ 7Xi xf ib ^ (Xic)& are 



100%. Masses are in GeV and the cross sections are in pb. The third row is the baseline point. 
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FIG. 8. The five regions in the supersymmetry parameter M2 — /i plane where different mass hierarchies occur. The three 
lightest neutralinos are denoted Xi,2,3i respectively. The dashed line is the locus of points scanned for the limits and is given 
by M2 = 0.89|/i| + 39 GeV. The dot is the baseline model described in the text. 
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FIG. 9. The limits on the cross section times branching ratio for SUSY production of 7&;Et events in the X2 
All production processes have been included; the dominant mode is the production of squarks and gluinos which decay to 
charginos and neutralinos. The overall branching ratio to the ^bfJ^ topology is approximately 30%. 
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FIG. 10. The limits on the X2X2 cross section in the X2 — > 7Xi SUSY model. 
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B. Gauge- mediated Model 



This is the second SUSY model |19| which can give substantial production of the signature jb^f. In this model the 
difference between the mass of the standard model particles and their SUSY partners is mediated by gauge (the usual 
electromagnetic, weak, and strong) interactions instead of gravitational interactions as in SUGRA models. SUSY 
is assumed broken in a hidden sector. Messenger particles gain mass through renormalization loop diagrams which 
include the hidden sector. SUSY particles gain their masses through loops which include the messenger particles. 

This concept has the consequence that the strongly-interacting squarks and gluinos are heavy and the right-handed 
sleptons are at the same mass scale as the lighter gauginos. A second major consequence is that the gravitino is very 
light (eV scale) and becomes the LSP. The source of b quarks is no longer the third generation squarks, but the decays 
of the lightest Higgs boson. If the lightest neutralino is mostly Higgsino, the decay xl hG can compete with the 
decays x? — > ZG and Xi — * jG- The Higgs decays to bb as usual. Since SUSY particles are produced in pairs, each 
event will contain two cascades of decays down to two x?'s, each of which in turn will decay by one of these modes. If 
one decays to a Higgs and one decays to a photon, the event will have the signature of a photon, at least one 6-quark 
jet, and missing Et- 

We will use a minimal gauge-mediated model with one exception. This MGMSB model has five parameters, with 
the following values: 

• A = 61 - 90 TeV, the effective SUSY-breaking scale; 

• M/A — 3, where M is the messenger scale; 

• N = 2 the number of messenger multiplets; 

• tan /3 = 3; 

• the sign of /i < 0. 

We wiU compute the MGMSB model using the GMSB option of ISAJET Q. We then re-enter the model using the 
MSSM options so that we can make one change: we set /i = — 0.75Mi. This makes the lightest neutralino a Higgsino 
so the branching ratio for Xi ~^ hG will be competitive with Xi — iG. We produce all combinations of xf' and Xji 
which are the only significant cross sections. We vary A which varies the overall mass scale of the supersymmetric 
particles. 



The model masses and branching ratios are given in Table XIV . The branching ratio is defined as the number of 
events with Xi iG and Xi ~^ hG divided by the number of events produced from all sources predicted by the 
model. We are using ISAJET with the CTEQ4L parton distribution function. The first point appears to have an 
unusually large efficiency because of other sources for b quarks which are not reflected in the definition of the signal 
branching ratio. We use the systematic uncertainties evaluated using the direct production of the X2 ~^ 7Xi model. 

Taking the two events observed, and convoluting with a 20% systematic uncert ainty gives an upper limit of 6.4 
events observed at 95% C.L. The final limits on this model are presented in Table |XV| and displayed in Figure |l|. 
Again, one can see that the experimental sensitivity is not adequate to set a mass limit (this time on the Xi mass) 
by several orders of magnitude. 
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M-± 

X, 


BR{xl ^ 7G) 


BR{x\ ^ hG) 


113 


141 


130 


90 


2 


132 


157 


147 


62 


18 


156 


178 


170 


33 


40 


174 


194 


186 


22 


50 



TABLE XIV. The models used in the hmits on the GMSB scenario. The Ughtest Higgs boson is 100 GeV. The masses are 
in GeV and the branching ratios are in %. 



At 


BR 


ffth X BR 


f95% lira X BR 


27 A 


3 


0.010 


0.27 


7.5 


20 


0.0402 


1.00 


8.4 


23 


0.0230 


0.89 


11.4 


18 


0.0111 


0.66 



TABLE XV. Efficiencies and Umits on direct production of Xf^Xj in the GMSB scenario. Branching ratios are not included 
in these efficiencies. The first row has an inflated efficiency due to the definition of the branching ratio. The units of Ae and 
the branching ratio are % and the cross sections are in pb. 
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FIG. 11. The limits on the cross section times branching ratio for SUSY production of 'ybfi^ events in the GMSB model. All 
production processes have been included. 



VI. MODEL-INDEPENDENT LIMITS 



As described in the introduction, there are several advantages to presenting limits of searches in a model-independent 
form. In the previous sections we derived limits on models of supersymmetry and presented the results as a limit 
on a cross section times branching ratio for a specific model, {a x BRY""' = N^""' /CAe, where is the 95% 

confidence level on the number of events of anomalous production and C is the integrated luminosity. We make a 
distinction between the acceptance. A, which is defined as probability that an object passes ry, and cuts, 
and the efficiency, e, which is the probability of events surviving all other sources of inefficiencies, such as photon 
identification cuts or 6-tagging requirements, which is detector-specific. The acceptance may be calculated from 
kinematic and geometric criteria alone, so an experienced worker in the field can compute it using only a Monte 
Carlo event generator program, while the efficiency requires access to the full detector simulation and, typically, 
multiple control samples. In our formulation, Num includes the degradation in sensitivity due to uncertainties on Ae, 
luminosity, and background subtractions, when they are included, as well as the statistical upper limit on the number 
of events. 

In the case of model-independent limits, there is no model to determine the efficiency and therefore we report 
a limit on (ct x BR x Aef"^™ = N^"^"^ / C. These limits, which are presented in the next section, do not have an 
immediate interpretation. (They do imply, however, a cross section range that we are not sensitive to, even with 
perfect efficiency.) In order to determine the meaning of these limits, in particular if a model is excluded or not, there 
must still be a mechanism for an interested physicist to calculate Ae for the model, and we develop three methods in 
the Appendix. 



A. Model— independent Limits on 'yhX Signatures 

The limit on {a x BR x Ae)'*™ = N^"^™- / C is described by reporting that two events are observed with an isolated 
photon with Et >25 GeV and \'q\ < 1, a SVX-tagged &-quark jet with Et >30 GeV and \'q\ < 2, A0(7- < 168° and 
> 40 GeV. The cuts are fully described in Section |ll[ The integrated luminosity for this sample is (85 ± 3) pb~^. 

The resulting 95% confidence level limit on {a x BR x Ae)'*™ for the ^hpi signature is then 0.069 pb. Adding the 
4% luminosity uncertainty we find the cross section limit increases to 0.070 pb. If we also add the 22% uncertainty in 
Ae from the WW limits (a typical uncertainty on an efficiency for this signature) discussed in the Appendix, we find 
the cross section limit increases 10% to 0.077 pb. This is the final model-independent limit on the signature ^hf]^. 
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The limit on the 76 signature before any JS^ requirement is 5.9 pb and the Hmit from the ^hJS^ signature from the 
98-event sample with > 20 GeV is 0.99 pb. 

The search for other objects in these events is described in Sections [IE and [VC, When we find no events, we can 
set a 95% confidence level limit on {axBRx eAY'"" of 0.038 pb assuming 4% uncertainty in the luminosity and 22% 
uncertainty in the efhcicncy This would apply to the searches for events with an additional photon, a muon or tau. 
For events with an additional electron, we observe one event and our limit becomes 0.057 pb. 

For events with N or more jets as shown in Figure we find the limits listed in Table XVl. 



iV 


(pb), > GeV 


4^° (pb), > 20 GeV 


1 


5.9 


0.99 


2 


2.5 


0.77 


3 


1.5 


0.50 


4 


0.52 


0.18 


5 


0.24 


0.10 


6 


0.14 


0.062 


7 


0.083 


0.038 


8 


0.062 


0.038 



TABLE XVI. The 95% confidence level limits on (a x BR x eA)'*™ in pb for events with TV or more jets, including the b jet. 



VII. CONCLUSIONS 



We have searched in 85 pb~^ of CDF data for anomalous production of events with a high-i?t photon and a 6-tagged 
jet. We find 1175 events with a photon with Et > 25 GeV and a 6-tagged jet with Et > 30 GeV, versus 1040±186 
expected from standard model backgrounds. Further requiring missing transverse energy > 40 GeV, in a direction 
not back-to-back with the photon (At/) < 168°), we observed 2 events versus 2.3± 1.6 events expected. In addition 
we search in subsamples of these events for electrons, muons, and tau-leptons, additional photons, jets and 6-quark 
jets. We conclude that the data are consistent with standard model expectations. 

We present limits on three current models of supersymmetry. The first is indirect production of chargino-neutralino 
pairs through squark and gluino production, where the photon is produced in ^2 ~^ 7Xi ^^'^ the 6-quark comes from 
the chargino decay into the light stop squark ib. A choice of favorable values of the parameters allows setting a 

lower mass limit on the gluino mass of 250 GeV. The second model is similar, but we look only at direct production 
of the X2xt pair- A cross section limit of ^ 7-10 pb is set, but is above the predictions for all X2 niasses so that no 
mass limit can be set. Lastly, a GMSB model is considered in which the photon comes from the decay x° — > 7G. 
Limits in the range 0.3-1.0 pb are set versus the mass of the xl, but again no mass limit can be set as the cross section 
predictions are lower than the limit. 

Finally, we present a model-independent limit of 0.077 pb on the production of events containing the signature 76 ^E^, 
and propose new methods for applying model-independent limits to models that predict similar broad signatures. 
We conclude that an experienced model-builder can evaluate whether model-independent limits apply to a particular 
model with an uncertainty of approximately 30%. 
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IX. APPENDIX: APPLICATION OF MODEL-INDEPENDENT LIMITS 



In the body of this paper, we present the hmits on specific models of new physics that predict the 76 signature, 
then rigorously calculate Ae for that model by using a Monte Carlo with a full detector simulation. We present our 
limits on (ct X Si?)"" = N^'"^/£Ae, or a parameter of the model such as the mass of a supersymmetric particle. 

A new paradigm, the signature-based or, equivalently, model-independent search may be an effective method for 
reporting the results of searches in the future. In this case, a signature, such as the photon and 6-quark jet addressed 
in this paper, is the focus of the search rather than the predictions of a particular model. 

There are several advantages to this approach 0J|,^ . 

1. The results are not dated by our current theoretical understanding. 

2. No a priori judgement is necessary to determine what is an interesting model. 

3. The results more closely represent the experimental observations and results will be presented in a form that 
can be applied to a broad range of models including those not yet imagined. 

4. The number of signatures is more reasonably limited than the number of models and model parameters. 

5. Concentrating on a particular model can tend to focus the search very narrowly, ignoring variations on the 
signature which may be just as likely to yield a discovery. 

6. Time spent on studying models can be diverted to systematically searching additional signatures. 

In order to reflect the data results more generally, in the body of this paper we also present a limit on (u x BR x 
^g-jhm _ j^jiirn I £ ^j^g signatures with no calculation of Ae. With limits presented this way, the collaboration itself, 
model-builders and other interested workers are no longer given limits on the physics models directly but now must 
derive the limits themselves. This has the potential for a wider application of the limits. In a pratical sense, it means 
the interested workers must calculate Ae for the model under study. 

In this Appendix, we present three methods to calculate Ae. These results together with the model-independent 
limits can be used to set limits on most models that predict events with the ")bp^ signature. 

The three methods are referred to as "object efficiencies" , "standard model calibration process," and "public Monte 
Carlo" . In the sections below we describe each in turn. In the following sections, we test these methods by comparing 
the results of each Ae calculation to the rigorously-derived Ae for the specific supersymmetry models. 



A. Object Efficiencies 



The first method for deriving Ae to use in conjunction with the model-independent limits is object efficiencies. 
The person investigating a model would run a Monte Carlo generator and place the acceptance cuts on the output 
which will determine the acceptance, A. The next step would be to apply efhciencies (simple scale factors) for the 
identification of each object in the signature, such as the photon or the 6-quark tag. This has the advantage of being 
being very straightforward and the disadvantage that correlations between the objects in the event are not accounted 
for. For example, a model with many jets would tend to have a lower efficiency for the photon isolation requirement 
than a model with few jets and this effect would not be reflected in this estimate of the efficiency. 

Using a sample oi Z ^ e+e~ events to measure the efficiencies of the global event cuts, we find the z < 60 cm cut 
is 92% efficient. The probability of finding no energy out-of-time is 98%. In this case the total global efficiency would 
be the product of these two efficiencies. In the discussion below, the efficiency of the identification of each object is 
often listed as efficiencies of several separate steps which should be multiplied to find the total efficiency. 

We can also use Z e+e~ events to measure the efficiency of the photon identification cuts. One electron from 
the Z is required to fire the trigger, but the second electron is unbiased with respect to the trigger. In addition the 
Z peak indicates the number of true physics events, ideal for measuring efficiencies. Which Z electron is required to 
pass which set of cuts (trigger or offline) must be effectively randomized to avoid correlations between the two sets 
of cuts. Requiring the cluster to be far from the boundary of the active area in the calorimeters is 73% efficient 
The trigger is 91% efficient, the identification cuts are 86% and the isolation requirement is 77% efficient. 

For the 6-quark efficiency we use a 70% probability that the jets from the event are contained in the SVX. (This 
would be 64% if the global event vertex was not already required to have z < 60 cm.) We add a 90% probability that 
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the jet was taggable (containing two reconstructed tracks in the SVX, passing pt cuts) and apply the pubHshed ]3^ ] 
tagging probabihty as a function of the jet Et which can be summarized as 



for Et < 18 GeV 
0.35 + 0.00277 * Et for 18 < i;* < 72 GeV 
0.6 for Et > 72 GeV 

The missing Et is found as the vector sum of the noninteracting particles in the event. As long as the missing Et 
is large, the resolution on the should not greatly effect the efficiency. 



In Section IV C, we searched the events in the 76 sample for additional leptons; here we present approximate object 
efficiencies for those cuts. These requirements and their efficiencies are borrowed from top-quark analyses ||ll|;|l2| as a 
representative selection for high-Et leptons. The efficiencies quoted here are measured in those contexts and therefore 
they are approximations in a search for new physics. In particular, the isolation efficiency is likely to be dependent 
on the production model. For example, if a model of new physics contained no jets, then the isolation efficiency is 
likely to be greater than that measured in top-quark events which contain several jets on average. 

For the electron search we require Et > 25 GeV and \r]\ < 1.0. Given that an electron, as reported in the output of 
the Monte Carlo generator, passes these acceptance cuts, the probability that the electron strikes the calorimeter well 
away from any uninstrumented region is 87%, the probability to pass identification cuts is 80%, and to pass isolation 
cuts is approximately 87% pl[ |. 

For muons we require pt > 25 GeV and |?7| <1.0. Given that the muon, as reported in the output of the Monte 
Carlo generator, passes these cuts, the fiducial acceptance of the muon detectors is 48%. Once the muon is accepted, 
the probability to pass identification cuts is 91%, and to pass isolation cuts is approximately 81%. 

Tau leptons are identified only in their one- and three-prong hadronic decays which have a branching ratio of 65%. 
(Tau semileptonic decays can contribute to the electron and muon searches.) We require that the calorimeter cluster 
has Et > 25 GeV and \t]\ < 1.2 and the object is not consistent with an electron or muon. Given that the r decays 
to a one- or three-prong hadronic decay mode and passes the Et and rj requirements, the probability that the tau 
passes identification and isolation cuts is approximately 57%. 



In Section IX D we apply these object efficiencies to the supersymmetry models and compare to the results of the 



rigorously-derived efficiency to test the accuracy of the results. 



B. Standard Model Calibration Process 



The second method for determining Ae for a model is the Standard Model process or efficiency model. In this method 
we select a simple physics model that produces the signature. The model is purely for the purpose of transmitting 
information about Ae so it does not have any connection to a model of new physics. Since it may be considered a 
calibration model, it does not have to be tuned and will not become dated. The interested model builder runs a 
Monte Carlo of the new physics and places acceptance cuts on the output, determining A, the same as the first step 
in the object efficiencies method. This result is then multiplied by the value of e which is taken to be the same as the 
value of e which we report here for the standard model process. 

We have adopted WW production as our efhciency model. One W is required to decay to ei' and we replace the 
electron with a photon before the detector simulation. The second W is forced to decay to bu, so the combination 
yields the signature jb^t- Since some efficiencies may be dependent on the Et of the objects in the event, we will 
vary the 'W mass to present this effect. A model-builder would then choose the efficiency that most closely matches 
the mass scale of the new physics models. 
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The Ae for this model is found using the same methods as used for the models of supersymmetry. From the 
difference in the observed efficiencies in Monte Carlo and data Z ^ ee samples, we use a 14% uncertainty on the 
efficiency of the photon ID and isolation. We use 9% for the 6-tagging efficiency uncertainty. The parton distribution 
function wc use is CTEQ4L. Comparing this efficiency to those obtained with the MRSDO' and GRV-94LO parton 
distribution functions, we find a standard deviation of 5%. Turning off initial- and final-state radiation increases 
the efficiency by 2% and 23% respectively and we take half these numbers as the systematic uncertainty. Varying 
the jet energy scale by 10% causes the efficienc y to c hange by 6%. We use an 4% uncertainty for the luminosity. In 
quadrature, the total systematic is 22%. Table XVI] summarizes the results. 



Mw 


75 GeV 


100 GeV 


125 GeV 


150 GeV 


175 GeV 


eA{%) 


0.85 


2.56 


4.86 


6.98 


8.12 


e{%) 


11.8 


10.7 


13.9 


15.4 


16.5 



TABLE XVII. Summary of the efficiencies found for the values of W mass used in the WW calibration model versus the 
value of the 'W mass. 
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In Section [XD we apply this method of calculating Ae to the supersymmetry models and compare to the results 
of the rigorously-derived efficiency to test the accuracy of the results. 



C. A Public Monte Carlo 



A Monte Carlo event generator followed by a detector simulation is the usual method for determining the efficiency 
of a model of new physics. However, there is usually considerable detailed knowledge required to run the simulation 
programs correctly so it is not practical to allow any interested person access to it. But if the simulation is greatly 
simplified while still approximating the full program, it could become usable for any worker in the field. The model- 
builder then only has to run this simple Monte Carlo to determine Ae. 

An example of this kind of detector simulation, called SHW, was developed for the Fermilab Run II SUSY/Higgs 
Workshop |33|. Generated particles are traced to a calorimeter and energy deposited according to a simple fractional 
acceptance and Gaussian resolution. A list of tracks is also created according to a simple efficiency and resolution 
model, and similarly for muon identification. The calorimeter energy is clustered to find jets. Electromagnetic clusters, 
together with requirements on isolation and tracking, form electron and photon objects. The tagging of 6-quarks is 
done with a simple, parameterized efficiency. At points where object identification efficiencies would occur, such as 
a cut on an electron shower profile, the appropriate number of candidates are rejected to create the inefficiency. 
The result is a simple list of objects that are reconstructed for each event. This method of determining efficiencies 
addresses the largest concern not addressed in the previous methods - the correlation of the characteristics of jets in 
the model with isolation requirements. We note that a highly parameterized Monte Carlo has obvious limitations. 

We have used the SHW program to compute efficiencies for the three models considered above. Since the program 
is tuned to provide the approximate efficiency of the Run II detector, we made a few minor changes to reflect the Run 
I detector. In particular, we changed the photon fiducial efficiency from 85% to 73% and the offiine efficiency from 
85% to 60%. We reduced the SVX acceptance along the z axis from 60 cm to 31 cm. Finally, we removed soft lepton 
6-tagging and added a 90% efficiency for the global event cuts. 

In the next section we use the public Monte Carlo to calculate Ae for the supersymmetry models and compare the 
results to the rigorously-derived efficiency to test the accuracy of the public Monte Carlo. 



D. Tests of the Model— independent Efficiency Methods 



In the body of this paper we have provided rigorous limits on several variations of three supersymmetry physics 
models that produce the signature of jb^^. In this section we apply the model- independent efficiency methods to 
the supersymmetry models. We can then compare the results with the rigorous limits to evaluate how effective it is 
to apply the model-independent limits to real physics models. 

In most cases we need to distinguish between acceptance and efficiency. Acceptance, indicated by A, we define as 
the probability for generated Monte Carlo objects to pass all geometric and Et cuts. For the 76^^ signature, with 
^( defined as the vector sum of neutrinos and LSP's, the cuts defining the acceptance of the signature are listed in 



Table XVIII 



7 


Et > 25 GeV 


\v\ <l-0 


6-quark 


Et > 30 GeV 


\v\ <2.0 


Additional Signatures 


Ept(j/, LSP) 


> 40 GeYAHl - ^t) < 168° 




e 


Et > 25 GeV 


\v\ <l-0 




Et > 25 GeV 


hi <l-0 


T 


Et > 25 GeV 


\V\ <1.2 


2'"' 7 


Et > 25 GeV 


\v\ <l-0 


2"=* 6 


Et > 30 GeV 


hi <2.0 


Jets 


Et > 15 GeV 


\v\ <2.0 



TABLE XVIII. The list of requirements on the output of a Monte Carlo generator which define the acceptance of a signature, 
A. The requirements on the photon and 6-quark jet above the double line are common to all signatures in this paper. When 
missing Et is required, as in all the supersymmetry searches and the tests of model-independent methods, both fi^ >40 GeV 
and A(/!>(7 — < 168° are required. The requirement is removed and other requirements are added to make specific 
subsamples. 
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Table XIX and Figure |12| list the results. The columns marked R are the efficiency times acceptance done rigorously, 
divided by the the same found using each of the model-independent methods. The difference of this ratio from 1.0 
one is a measure of the accuracy of the approximate methods compared to the rigorous method. 
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18 


54.7 
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1.26 


1.22 


1.11 


q, g production 

Ms=Mg 


185 


30 


17.0 


1.97 


0.91 


0.68 


0.48 


210 


30 


22.0 


2.98 


1.04 


0.73 


0.90 


235 


30 


24.0 


3.23 


1.01 


0.68 


0.90 


260 


30 


24.5 


2.69 


0.82 


0.52 


0.75 


285 


30 


19.7 


2.16 


0.84 


0.48 


0.72 


-> 7X1 
q, g production 

Ms = M,± 

^2 


110 


100 


13.5 


0.93 


0.54 


0.54 


0.59 


130 


100 


12.6 


1.41 


0.88 


0.80 


0.87 


140 


100 


14.8 


1.29 


0.68 


0.60 


0.66 


150 


100 


13.7 


1.34 


0.77 


0.65 


0.78 


170 


100 


11.5 


1.27 


0.85 


0.68 


0.65 



TABLE XIX. The results of comparing the methods of calculating Ae using the model-independent methods and the 
rigorously-derived Ae. Each row is a variation of a model of supersymmetry as indicated by the label in the first column and 
the mass of a supersymmetric particle listed in column two (GeV). The column labeled A is the acceptance of the model in % 
and the next column is the rigorously-derived At. The columns labeled with R are the ratios of the rigorously-derived Ae to 
Ae found using the model-independent method indicated. 
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FIG. 12. The ratio of the efficiencies obtained wi th fu ll detector simulation to those obtained with the model-independent 
methods. The x axis is the row number from Table XIX 



E. Conclusions from tests 



There are several notable effects apparent immediately from Table XIX and Figure |lj. The first is that the 
comparison of efficiencies for one model point fares especially poorly. This occurs when the branching ratio for the 
model is very small (2%). When the events do not contain many real photons and b quarks, the small number of 
objects misidentified as photons and 6-quarks becomes important. For example, jets may pass photon cuts and c 
quarks may be 6-tagged. When this occurs, the full simulation will be more efficient than a method which specifically 
requires that the Monte Carlo generate an isolated photon or 6-quark in order to accept the event. This is true of the 
object efficiencies method and the efficiency model method. We note that the public Monte Carlo method does allow 
misidentification and so it does not show this large mismatch. We can conclude that when the branching ratios are 
small, the public Monte Carlo method is vastly superior to the others. 

In the object efficiency method, the acceptance of the signature is computed by running the Monte Carlo without 
a detector simulation. As each object in the signature is identified and passes acceptance cuts, the individual object 
efficiencies are multiplied. These object efficiencies which may or may not be Et- or ry-dependent, are listed in Section 



IX A. In this test, these efficiencies are typically well-matched to the rigorously-derived efficiencies. The average of 
Robj over all models except the first, is 0.88 ± 0.21 where 0.21 is the RMS computed with respect to 1.0, the ideal 
result. 

In the efficiency model method, we generate a Monte Carlo model that is not related to a search for new physics 
but produces the signature of interest. For the signature of jb^^., we generated WW {'^v){bu). The efficiency 
model results are also optimistic, the average is a ratio of 0.74 ±0.35 where again the uncertainty is actually the RMS 
with 1.0, the ideal result. We found that the difficulty of applying this method was in choosing the mass scale. For 
example, we chose to match the 'W mass to the mass in the direct production of the X2 ~* 7Xi model. However, 
the photon comes from a secondary decay and the effect of the LSP mass compared to the massless neutrino causes 
the Et of the 7 and b to be poorly matched to the Et of these objects in the WW model. 

In the public Monte Carlo method, we compute the efficiency using SHW, a highly-parameterized, self-contained 
Monte Carlo. In general, results here are somewhat optimistic with the average ratio to the rigorous total efficiency 
being 0.77 ± 0.28, where the uncertainty is the RMS computed with respect to 1.0, the ideal result. 

For completeness we also include the ratio of the simple acceptance to the rigorous acceptance times efficiency. The 
average ratio is 0.12 ± 0.87. 

The methods for calculating efficiency without access to the full detector simulation are accurate to approximately 
30% overall. They tend to underestimate Ae by 10-25% but the result for individual comparisons varies greatly. These 
uncertainties are larger than, but not greatly larger than, a typical uncertainty in a rigorously-derived efficiency, which 
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is 20%. 

We conclude that to determine if a model is easily excluded or far from being excluded by the data, the model- 
independent methods are sufficient. If the model is within 30% of exclusion, no conclusion can be drawn and the 
efficiency should be rigorously-derived. 
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